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ABSTRACT 



We analyze the long-term evolution of the fluxes of six active galactic nuclei (AGN) - 0923-1-392, 3C 1 1 1 , 3C 273, 3C 345, 3C 454.3, 
and 3C 84 ~ in the frequency range 80-267 GHz using archival calibration data of the IRAM Plateau de Bure Interferometer. Our 
dataset spans a long timeline of » 14 years with 974 - 3027 flux measurements per source. We find strong (factors a:2"8) flux variability 
on timescales of years for all sources. The flux density distributions of five out of six sources show clear signatures of bi- or even 
multimodality. Our sources show mostly steep (a x 0.5 - 1), variable spectral indices that indicate outflow dominated emission; 
the variability is most probably due to optical depth variations. The power spectra globally correspond to red-noise spectra with five 
sources being located between the cases of white and flicker noise and one source (3C 111) being closer to the case of random walk 
noise. For three sources the low-frequency ends of their power spectra appear to be upscaled in spectral power by factors a; 2-3 with 
respect to the overall powerlaws. In two sources, 3C 454.3 and 3C 84, the 1.3-mm emission preceeds the 3-mm emission by a;55 
and a;300 days, respectively, probably due to (combinations of) optical depth and emission region geometry effects. We conclude that 
the source emission cannot be described by uniform stochastic emission processes; instead, a distinction of "quiescent" and (maybe 
multiple) "flare" states of the source emission appears to be necessary. 

Key words. Galaxies: active — Quasars: general — Radiation mechanisms: non-thermal 



1. Introduction 

Active galactic nuclei (AGN) have been studied extensively in 
the wavelength range from cm-radio to y-rays in the last decades 
(see, e.g., Kembhavi & Narlikar [T999l or Krolik [T999l and refer- 
ences therein for a review). There is overwhelming observational 
evidence that their emission originates from accretion onto su- 
permassive black hole s (SM BH) with masses M, ^ IO^-^Mq 
[ (e.g., Ferrarese & Ford l2005 l and references therein). 

The viewing-angle unification standard model of AGN 
(Lawrence' 1987'; see also Urry & Padovani 1995) has long been 
established, although recent studies emphasize the impact of 
' galaxy mergers on AGN fueling and luminosity (e.g. Hopkins 
' et al. 120061 ). However, the details of AGN activity are still not 
well understood. One example is the nuclear emission; propos- 
als for its origin range from shocks in continuous (e.g., Marscher 
& Gear . 1985 ) or discontinuous (e.g., Spada et al. 12001!) jets to 
orbiting plasma "hotspots" (e.g., Abramowicz et al. I1991I I or 
plasma density waves (e.g., Kato'2000') in accretion disks. 

One possibility to constrain the AGN physics is offered by 
the analysis of characteristic timescales and the noise properties 
of the emission. Several studies aimed at identifying character- 
istic timescales, including possible quasi-periodic oscillations. 



from tens of minutes (using X-ray observations; e.g. Benlloch 
et al. 12001! and using mm-radio observations; e.g. Schodel et 
al. !2007! l to tens of years (based on cm- and mm-radio moni- 
toring data; Hovatta et al.'2 007!l2008! l. Other works focused on 
the power spectral statistics of the emission which was found 
to globally follow red-noise type powerlaws (e.g. Lawrence & 
Papadakis [T993l Uttley et aL l2n02l l. 



* This study is based on observations carried out with the IRAM 
Plateau de Bure Interferometer. IRAM is supported by INSU/CNRS 
(France), MPG (Germany), and IGN (Spain). 



The studies by Hovatta et al. (^SOOTl l2008l l used radio 
data in the frequency range 22-87 GHz spanning timelines of 
*25 years. A key result of their study was that AGN show 
large-amplitude activity on timescales >6 years. Therefore, 
analysing the long-term behaviour of AGN requires monitor- 
ing on timescales of >10 years. Building up on this conclusion, 
we have analyzed archival flux monitoring data from the IRAM 
Plateau de Bure Interferometer (PdBI) of six AGN - 0923-H392, 
3C 111, 3C 273, 3C 345, 3C 454.3, and 3C 84 - in the (obser- 
vatory frame) frequency range 80-267 GHz. Given the range of 
redshifts z ~ 0.02 - 0.86 of our sample, this corresponds to a 
range of rest-frame frequencies of ai82^97 GHz. Our dataset 
comprises several thousand measurements for each source, cov- 
ering timelines of a; 14 years. We note that our study comple- 
ments earlier AGN mm/radio surveys carried out with the IRAM 
Pico Veleta observatory (Steppe et al [19881 [19921 [19931 Reuter 
etaL [T997l i. 
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Table 1. Physical properties and observation journal for our six target AGN. J2000 coordinates, source types, and redshifts are taken 
from the NED, black hole masses M, are from the references given. The kpc-scale morphology has been taken from the MOJAVE 
database (Lister et al. 2009). We also give the total monitoring time T and the number of measurements A^ for the 1.3-mm (here 
labeled 1-mm for brevity), 2-mm, and 3-mm bands, respectively. 



Object 



R.A. 



Dec 



Type" Redshift M.IW^Mq] kpc-morpliology T [yr] A^i, 



A'2, 



A'3„- 



0923+392 


09:27:03 


+39:02:21 


FSRQ 


0.695 


n'- 


1 -sided halo 


13.7 


574 


42 


1019 


3C 111 


04:18:21 


+38:01:36 


Syl 


0.049 


36'^ 


2-sided jet 


12.9 


398 


7 


569 


3C273 


12:29:06 


+02:03:09 


FSRQ 


0.158 


8" 


1 -sided jet 


13.8 


1053 


81 


1893 


3C345 


16:42:58 


+39:48:37 


FSRQ 


0.593 


18'' 


1 -sided halo 


13.9 


587 


57 


1099 


3C 454.3 


22:53:57 


+ 16:08:54 


FSRQ 


0.859 


7" 


1 -sided jet 


13.8 


835 


116 


1677 


3C84 


03:19:48 


+41:30:42 


Sy2 


0.018 


19c 


2-sided halo 


13.9 


501 


56 


990 



" "Sy" = Seyfert galaxy, "FSRQ" 

*> Liu etal. (2006) 

' Marchesini et al. J2004t 



: flat spectrum radio quasar 



2. Observations and Data Processing 

The PdBI is able to observe in three atmospheric windows lo- 
cated around wavelengths of 1.3 mm, 2 mm, and 3 mm. Each 
of these bands covers a continuous range of frequencies that are 
available for observations; these ranges are 201-267 GHz for 
the 1.3 mm band, 129-174 GHz for the 2 mm ba nd, and 80- 
1 16 GHz for the 3 mm band (Winters & Neri|20T3). 

Until the end of 2006, observations were made using dual- 
band, single-polarization receivers. One frequency in the 1.3- 
mm band and one frequency in the 3-mm band could be se- 
lected for simultaneous observations; the 2-mm band was not 
reachable. At the beginning of 2007, new single-band, dual- 
polarization receivers were installed. Since then it is possible to 
observe a single frequency located in a given band simultane- 
ously in both linear polarizations. At the beginning of 2008, the 
2-mm band became available for observations. 

Flux monitoring data for AGN are available from October 
1996 on as a by-product of regular observatory operations. The 
PdBI uses the most luminous (Sy ^ 3 Jy in the 3-mm band) ac- 
tive galactic nuclei as calibrators for the spectral response of the 
instrument. The sources we present here are the most luminous 
targets regularly available (especially in terms of declination) as 
spectral calibrators. Therefore at least one of these sources has 
been observed in almost any individual observing project, with 
each project lasting typically few hours. This leads to a dense 
sampling over the entire observation timeline of > 1 3 years; gaps 
in the lightcurves are due to technical works at the interferometer 
and seasonal sunavoidance constraints. For our final sample of 
sources to be analyzed we picked all objects for which we have 
at least 1000 measurements over the entire observing timeline 
without major (>1 year) interruptions of the monitoring. These 
were 0923+392, 3C 1 1 1, 3C 273, 3C 345, 3C 454.3, and 3C 84. 

Antenna temperatures are converted into physical flux den- 
sities by using empirical antenna efficiencies as conversion fac- 
tors. Those factors are functions of frequencies and are located 
in the range from x!22 Jy/K (for the 3-mm band) to =s37 Jy/K 
(for the 1.3-mm band). In order to estimate the systematic un- 
certainties of our dataset, we analyzed observations of the radio 
continuum star MWC 349A (e.g. Tafoya et al. 2004) that were 
obtained and calibrated like the AGN measurements. Due to its 
well-known flux properties MWC 349 serves as a "standard can- 
dle" for the PdBI (Krips et al. in prep.). From the systematic scat- 
ter of the MWC 349 lightcurves we conclude that the systematic 
relative uncertainties of our AGN observations are a!lO% for the 
2-mm and 3-mm bands and a: 16% for the 1 .3-mm band. We note 



that the systematic accuracy of our data is limited by the auto- 
matic monitoring and calibration; when calibrating observations 
individually and interactively while using MWC 349 as refer- 
ence, systematic errors can be reduced to <5% (Krips et al. in 
prep.). 

Before analysis, we had to remove erroneous values from 
the flux data. Errors are caused by either of two effects: (1) poor 
signal-to-noise ratio (S/N), or (2) failure of the flux calibration. 
We caught effect (1) by rejecting data with statistical flux er- 
rors larger than 0.1 Jy or S/N<10. Statistical errors of individ- 
ual measurements used for analysis are thus <0. 1 Jy. Effect (2) 
is more difficult to treat as it leads to a systematic deviation of 
the observed from the true flux values. Especially, it leads to a 
clustering of data points in the range 0-1.5 Jy (details depend- 
ing on the source), well separate from, and systematically be- 
low, the long-term lightcurves. In order to remove this effect, we 
rejected flux measurements with fluxes located below the "di- 
vides" that separate the long-term lightcurves from the clusters 
of erroneous points. Additionally, we checked the statistical be- 
haviour of the data. For this, we analyzed the differences be- 
tween the lightcurves and smoothed versions of the lightcurves, 
using smoothing by running medians. We found that the distribu- 
tions of the residuals usually follow Gaussian distributions with 
extended tails. We regarded the data points located in the tails 
as "suspicious" outliers and rejected all data located more than 
2.5cr away from the centers of the distributions. This conserva- 
tive approach rejected between 2% and 16% of the data of a 
given lightcurve. 

The data processing as well as the subsequent analysis made 
use of the IRAM Grenoble data processing software GILDAq3 
and the MPE Garching data processing software DPUSEPQ. 

After processing the data, we were left with large numbers 
of measurements ranging from 974 to 3027 per target, corre- 
sponding to one measurements every 1.5^.3 days. Our flux 
data base covers a timeline from September 1996 (for the 2- 
mm band: February 2008) to August 2010. In Table [1] we pro- 
vide an overview over our target AGN and individual observa- 
tion statistics. We collected the source types and redshifts from 
the NASA/IPAC Extragalactic Database (NED) and the kilo- 
parsec scale morphology from the MOJAVE database (Lister et 
al. 120091 ). 



' Developed and maintained by the IRAM GILDAS team, see 
http://www.iram.fr/IRAMFR/GILDAS/ 

- Developed and maintained by Thomas Ott (MPE Garching), see 
http://www.mpe.mpg.de/~ott/dpuser/index.html 
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3. Analysis 

3.1. Lightcurves 

In Fig. [T]we present the lightcurves of our six target AGN. For 
each source we give flux densities 5 v vs. time t separately for the 
1.3-mm, 2-mm, and 3-mm bands as far as available. 

In order to quantify the variability timescales, we approxi- 
mated the most prominent - referring to smooth, continuously 
sampled variations in flux by factors >2 - variations of each 
lightcurve as exponential rises or decays like 



Syoce 



±1/1. 



(1) 



with tf, being the evolution timescale and the sign of the exponent 
being positive for exponential rises and negative for decays. We 
attempted those approximations only for parts of lightcurves that 
were sufficiently smooth for being described by a single analyt- 
ical function and well sampled. These limitations prevented us 
from applying this approach to the lightcurve of 3C 111 which 
shows a combination of irregular sampling and fast variability. 
For the other sources, we give the results in the corresponding 
diagrams of Fig. [T] 

3.2. Flux Distributions 

In Fig. |2] and [3] we present for each source the flux distributions 
for the 1.3-mm and 3-mm bands. The 2-mm lightcurves do not 
contain enough data for calculating corresponding distributions. 
In order to prevent distortions introduced by irregular sampling, 
we binned our data in time. The size of the time bins for a spe- 
cific lightcurve is given by 



At = 2T/N 



(2) 



Here T is the total observation period, A^ is the number of data 
points. For the special case of regular sampling. At corresponds 
to the inverse of the Nyquist frequency. 

As the resulting histograms show somewhat complicated 
morphologies, we investigated the minimum number of smooth 
components necessary for a reasonable description of their pro- 
files. For this we approximated each distribution as a superpo- 
sition of two (0923+392, 3C 111, 3C 345) or three (3C 273, 
3C 454.3, 3C 84) Gaussian distributions 



G(S)^aoexp[-{S -Sof/cr^] 



(3) 



with S denoting the flux. So being the center position, cio being 
the amplitude at So, and cr^ denoting the variance. The choice 
of Gaussian distributions was governed by the fact that a priori 
such distributions are most likely to occur in case of uniform 
stochastic emission processes. 

In three cases (0923H-392 at 1.3 mm, 3C 111 at 1.3 mm and 
3 mm) we also attempted to describe the flux histograms by log- 
normal distributions 



MS) 



1 



^fhio-S 



exp[-(ln5 -jifllo-^] 



(4) 



where In denotes the natural logarithm and <j and ji are the 
scale parameter and the location parameter, respectively. As log- 
normal distributions are characteristic for multiplicative stochas- 
tic processed they may be expected in AGN flux distributions 
(compare, e.g., Dodds-Eden et al. 1201 11 1. 



' This is a simple consequence of the central limit theorem which 
states that sums of random numbers are normally distributed. As the 



3.3. Spectral Indices 

As we observed our target AGN over a range of frequencies 
located in three spectral bands, we were able to analyze their 
mm-spectroscopic properties. This analysis was however com- 
plicated by fast intrinsic variability as well as by the fact that the 
spectral information partially needed to be derived from non- 
simultaneous observations. We therefore inserted all measure- 
ments within selected time windows into flux-frequency dia- 
grams and derived time-averaged spectral indices a via means 
of least-squares fits. Spectral indices were defined via the rela- 
tion 



(5) 



For each source we present the spectral index as function of 
time in Fig. |4] We selected the time windows such that a given 
window covers a phase of similar flux levels, meaning a "quiet" 
phase, an outburst, or a low-flux regime between two outbursts. 

In order to test possible coiTelations between spectral index 
and source flux, we present a as function of the 3-mm flux for 
each source in Fig.|5] For this diagram we use the same binning 
in time as for Fig.|4] In order to preserve some time information, 
we connect and label the data points according to their order in 
time. 



3.4. Power Spectra 

As the lightcurve of our sources show variability on many 
timescales, we performed a power spectrum analysis in order 
to quantify this variability. For each lightcurve we computed a 
normalized Scargle periodogram 



A,= 



1 

2^ 



{ZjSj cos Inftj) {ZjSj sin Inftj] 
ZjCOS^2j:ftj "^ Zj&in^2nftj 



(6) 



(Scargle Fl982l l from our data. Here Af is the periodogram am- 
plitude evaluated at frequency /, S j and tj are the flux value 
and the time of the j-th data point, respectively, and cr^ de- 
notes the variance of the data. The base frequencies are / = 

/min, 2/mm,3/;nin,..., /max with /min = 1/7^, /max = N/(2T),T 

being the total observation period, and A^ being the number 
of data points. For the special case of regular sampling, /max 
corresponds to the Nyquist frequency. The timelines and sam- 
pling of our observations (see Table [TJ lead to typical values of 
/min ~ 2 X 10"^ days"' and /^ax ~ 0.2 days"' . Before calculating 
a periodogram we subtracted the average source flux from the 
data in order to remove the zero-frequency power. In order to 
exclude systematic errors due to spectral eff'ects, we performed 
the periodogram analysis separately for the 1.3-mm and 3-mm 
bands of each source (for the 2-mm band the amount of data is 
insufficient for a comparison to the other bands). 

The power spectra we find all show a global decrease of 
spectral power with increasing frequency, i.e. they correspond 
to red-noise spectra. In order to quantify this behaviour, we de- 
rived powerlaw indices by fitting straight lines to the data in log- 
space. The spectral index /3 of the power spectra is defined via 



y 



f 



-p 



logarithm of a product of random numbers equals the sum of the log- 
arithms of the individual random numbers, it follows (under certain 
weak conditions) that the logarithms of products of random numbers 
are normally distributed (e.g., MontroU & Shlesinger |1982| and refer- 
ences therein). 



S. Trippe et al.: The long-term millimeter activity of active galactic nuclei 



In several cases we find that the low-frequency ends of the 
power spectra follow the average powerlaw slopes but are lo- 
cated systematically above the best-fitting powerlaw models. In 
those cases we approximated the affected spectral ranges as up- 
scaled (by factors ^2-3) versions of the overall power spectra. 
We present our results in Figs.|6]and|71 



by the range [max{DCF), max(DCF) - TxTmaxiDCF)], thus cover- 
ing all values of DCF that are in agreement with the maximum 
of the DCF on a 3cr confidence level. We probe the null hy- 
pothesis "the 1.3-mm and 3-mm lightcurves are simultaneous". 
Accordingly, we regard the lightcurves as simultaneous as long 
as DCF(t = 0) > max(DCF) - 'icTmaxi.DCF)- 



3.5. Time Offsets among Spectral Bands 

The long timeline of observations and the good time sampling of 
our flux data make it possible to test if the lightcurves of the var- 
ious spectral bands are simultaneous or if time lags are present. 
In order to quantify the presence or absence of time offsets, we 
computed for each source a discrete correlation function (DCF) 
according to the scheme proposed by Edelson & Krolik (I1988I I. 
For two discrete datasets {a,), {bj} with averages a, b one com- 
putes the unbinned discrete correlations 



UDCFij(Mj) = 



(fl; - a)(bj - b) 



(7) 



for all pairs {ai,bj) with time differences Af,j. The cr^^ are the 
variances, the 6a,h are the mean measurement errors of {a,), {bj}. 
Accordingly, the terms o-^i, - S^i, are the variances corrected 
for artificial broadening due to finite measurement errors; this is 
necessary in order to preserve the correct normalization in noisy 
data. The actual DCF for a time offset r is obtained via averag- 
ing over all A^ UDCFij{Atij) for which Af,j is located within a 
selected bin At (i.e., t - At/2 < Ar,y <t + At/2): 



t+At/2 



DCF{T)^- 2 UDCFijiAtij) 



(8) 



t-At/2 



By definition, DCF(t) is located in the range [-1, -Hi], with (-1) 
+ 1 indicating perfect (anti-)correlation and corresponding to 
the absence of any correlation. The position of the maximum of 
the DCF indicates the time offset between the lightcurves. The 
statistical uncertainty of DCF(t) is given by the standard error 
of mean 



o-dcf(t) - 



1 



N-l 



r+Ar/2 

2 {UDCFij-DCFiT)f 

t-At/2 



1/2 



(9) 



For each of our six target AGN, we examined ranges of r suf- 
ficient for covering the largest time offsets to be expected, either 
T ± 365 days or r ± 730 days. In the special case of 3C 454.3, 
we additionally probed two subsets of data with t + 180 days 
(see Sect. 4.5). In order to preserve high time resolution while 
excluding the risk of creating oversampling artefacts, we used 
sampling times At = Tin^m/Numi,, with Ti^,^ (A^imm) being the 
total observing time (number of measurements) obtained in the 
1.3-mm band which is for all sources the wavelength band with 
the sparser sampling (see Table[T]). We note that the effective res- 
olution depends only weakly on the choice of bin sizes: smaller 
At increase the formal resolution but also increase the errors as 
they cover smaller numbers of UDCFjj per bin (see Eq. |9ll, and 
vice versa. 

We present the results of our analysis in Fig. [8] In our con- 
vention, a positive (negative) time lag means that the 3-mm 
lightcurve preceeds (follows) the 1.3-mm lightcurve. In each di- 
agram, we indicate the ranges (as grey-shaded areas) covered 



4. Results 

4.1. 0923+392 (4C +39.25) 

The 3-mm (1.3-mm) flux densities of this flat spectrum radio 
quasar (FSRQ) dropped from x6 Jy (^3 Jy) in 1996 to x3 Jy 
(~1.5 Jy) in 2003 and have been rising again since then to 
present-day levels of ^4.5 Jy (~3 Jy). For the flux density de- 
crease from 1996 to 2003 we find characteristic (exponential) 
timescales of a;5-7 years. The 3-mm flux distribution shows 
a double peak profile that can be described by two Gaussian 
components separated by a; 1.6 Jy. The 1.3-mm flux distribution 
shows an asymmetry indicative of a second component but no 
significant secondary peak. We therefore attempted to describe 
the histogram by a log-normal distribution according to Eq. |4] 
However, we did not find a satisfactory solution, with the best 
fitting curve showing ax^/dof = 13.8. 

The spectral index a is a;0.8 in the first years of monitoring 
before 2000 and steepens to a; 1 around 2002. The index then flat- 
tens again until reaching a!0.6 in the time window 2007-2010. 
There appears to be no correlation between spectral index and 
3-mm flux. Even though we find the steepest spectral indices for 
the lowest source fluxes, we find similar variations in a (by aiO.2) 
within the same flux ranges. 

The 1.3-mm and 3-mm lightcurves are simultaneous within 
the 3cr confidence interval. The power spectra of the 3-mm and 
1.3-mm flux measurements can be well described as simple red- 
noise spectra with fairly flat powerlaw indices /3 ~ 0.4. Within 
the errors, the values for /? are in agreement for both frequency 
bands. 



4.2. 3C 111 

This Seyfert 1 galaxy shows episodic strong variability on 
timescales of years. The source seems to spend most of the time 
in a low-activity or "quiescent" state with fluxes in the range 
a; 1-3 Jy. In addition, we can identify at least five pronounced 
spikes in the Ughtcurves in 1997, 2001, 2006, 2008, and 2009. 
These outbursts reach flux densities of ai6-15 Jy, thus exceeding 
the quiescent level by factors 3-6. The outbursts appear to rise 
and decay on fast (<1 yr) timescales; unfortunately, the irregular 
sampling of the lightcurves prevents a reliable quantitative anal- 
ysis. The 3-mm (1.3-mm) flux distributions can be described by 
smooth, asymmetric profiles composed of two Gaussian compo- 
nents for fluxes <8 Jy (<5 Jy) plus sparsely populated tails rang- 
ing out to a;14 Jy (alO Jy). Except for the tails, it is also possible 
to describe the histograms with (unimodal) log-normal distribu- 
tions. In both cases we find best-fit solutions with x^ /do f a 1.6. 
Although this is a relatively high value, the log-normal distribu- 
tion appears to be an acceptable description of the data. Within 
errors, the 1.3-mm and 3-mm lightcurves are simultaneous. 

For most of the monitoring time the spectral index remains 
at a X 0.5 (within errors). The only exception appears around 
1997 when the spectral slope is a a 0.8; this phase coincides 
with the first outburst. In the flux-a diagram the various points 
do not indicate a trend. 
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The power spectrum of the 3-mm (1.3-mm) flux data glob- 
ally follows a red-noise powerlaw spectrum with /? a; 1 .0 (/? ~ 
0.6). The low-frequency end (/ ^ 0.002 - 0.01 days') of the 
3-mm power spectrum is located systematically above the best- 
fitting overall powerlaw; this region can be described as an up- 
scaled (by a factor »2) version of the global red-noise spectrum. 

4.3. 3C 273 

At the beginning of the monitoring, the 3-mm (1.3-mm) 
lightcurves of this FSRQ rose fast towards a strong outburst 
peaking in 1997 at levels of ^38 Jy (=^30 Jy). From then on, 
the flux densities steadily decreased down to minimum values 
of =s6 Jy (x4 Jy) in 2004 with an exponential decay timescale 
of ^3 yr (measured for the decrease in 1997-2000). Since about 

2006, 3C 273 has been in a state of high variability, with fluxes 
fluctuating in ranges ailO-27 Jy (»6-20 Jy). The flux distribu- 
tions follow complex, double-peaked profiles; they can be ap- 
proximated by superpositions of three Gaussian distributions 
peaking at a;8 Jy, -18 Jy, and a;28 Jy (~4 Jy, »10 Jy, and a;16 Jy). 
Within eiTors, the 1.3-mm and 3-mm lightcurves are simultane- 
ous. 

The spectral index varies substantially over the 14 years of 
observation, starting at c « 0.3 around 1997, steepening up to 
a K 0.8 in 2000-2003 and flattening again to o- a; 0.4 after 

2007. There is a global trend towards flatter indices at higher flux 
levels, though no clear correlation: on the one hand a remains at 
similar levels over wide flux ranges, on the other hand strong 
variations in a occur within the same flux ranges. 

The power spectra of the 3-mm and 1.3-mm Ughtcurves fol- 
low red-noise powerlaws with identical (within errors) indices 
y6 a 0.45. For frequencies / < 0.04 days"' the spectral power is 
located systematically above the best-fitting overall powerlaw; 
these ranges can be described as varieties of the overall power 
spectra upscaled by factors ^3 (3 mm) and ^2 (1.3 mm), respec- 
tively. 

4.4. 3C 345 

The 3-mm (1.3-mm) lightcurves of this FSRQ start at flux den- 
sities of ^4 Jy (=s2 Jy) in 1997 and reach their highest values of 
ss7 Jy (~4 Jy) in mid- 1998. From then the fluxes decrease over- 
all down to ^3 Jy (~2 Jy) in the time range 2004-2008. Around 
2009 an outburst occurs, reaching flux levels up to a!6.5 Jy 
(^4.5 Jy) and slowly decreasing since then. The 2009 outburst 
raises in 2008-2009 with an exponential timescale of a 1 yr. The 
3-mm (1.3-mm) flux distributions follow double-peaked profiles 
we approximated as superpositions of two Gaussians peaking at 
s»3 Jy and aS.S Jy (^2 Jy and a!3.5 Jy). Within errors, the 1.3- 
mm and 3-mm lightcurves are simultaneous. 

Within errors, the spectral index remains on very similar lev- 
els in the range a a 0.6 - 0.8 for most of the monitoring time. 
After 2007 - the time frame of the most recent outburst - the 
spectrum flattens somewhat towards a a 0.4. We note that this 
somewhat steep index at mm-wavelengths is consistent with the 
classification of 3C 345 as FSRQ: the classification as FSRQ 
is based on cm-observations, whereas at mm-wavelengths the 
spectrum bends towards a > 0. Ina5v-Q' diagram the data 
points move with time on an almost circular trajectory, leaving 
no room for a correlation between flux and spectral index. 

The power spectra can be described as simple red-noise spec- 
tra, with powerlaw indices/? ~ 0.5 for the 3-mm and/3 ~ 0.4 for 
the 1.3-mm lightcurve, respectively. 



4.5. 3C 454.3 

This flat spectrum radio quasar shows the most spectacular 
lightcurves among the sources of our sample. For the first nine 
years (1996-2005) of monitoring the activity of this source was 
restricted to variability on levels of a;2-8 Jy. This changed dra- 
matically in 2005 when the first of at least three major outbursts 
occured, peaking at flux densities of a!35^0 Jy. The series of 
adjacent outbursts has continued until the present day, making 
3C 454.3 one of the brightest extragalactic mm/radio sources 
cuiTently observed (see also, e.g., Jorstad et al. 1201 Oi l. All out- 
bursts show short exponential rise and decay times in the range 
~0. 1-0.5 years. The 3-mm (1.3-mm) flux distribution is dom- 
inated by a Gaussian peak centered at ^6 Jy (~4 Jy) with an 
extended tail ranging out to ^40 Jy (^35 Jy). 

The spectral index has been fluctuating around a a 0, with 
the most extreme values being a a -0.4 and a a 0.5, respec- 
tively. A somewhat steep value of a x 0.5 can be observed for 
the time 1996-2000. Notably, the most extreme change (from 
a X -0.4 to ff a 0.4) occurs at the time of the peak of the first 
bright outburst in 2005. The data almost fill the entire 5 y -a area 
observed, leaving the impression that flux and spectral index are 
uncorrected. 

The 3-mm (1.3-mm) power spectrum globally follows a red 
noise powerlaw with index /3 ~ 0.7 (J3 ~ 0.5). However, for 
frequencies / < 5 x 10"^ days"' the power spectrum seem to 
be an upscaled - by a factor a 3 (^2) - version of the average 
powerlaw. 

The discrete correlation function points to time delays t 
significantly different from zero in terms of the 3cr confidence 
threshold, located in the range t a -20... - 110 days. The neg- 
ative delay indicates that the 1.3-mm emission preceeds the 3- 
mm emission. However, as 3C 454.3 shows a dramatic transition 
from "quiescent" to "flaring" around epoch 2005.2, the curve 
presented in Fig. [8] (bottom left) might be misleading as the tim- 
ing of the emission might be different in the two states. We there- 
fore computed in addition one DCF for each of the two activity 
phases, one covering all data from 1996.8 to 2005. 1, one includ- 
ing all measurements obtained between 2005.2 and 2010.6. The 
results are summarized in Fig.|9] We note that for the "quiescent" 
phase (before 2005.2) the most probable time delay - on a low 
absolute level, DCF(Q) ^ 0.6 - is consistent with zero, whereas 
we find for the "flaring" phase (after 2005.2) r ^ -15... - 80 
days, with the maximum of the DCF located at t a -55 days. 
Again, the negative delay indicates that the 1.3-mm emission 
preceeds the 3-mm emission. 

Finally, we probed the impact of the spectral time delay on 
the calculation of the spectral index a. For this, we repeated 
the calculation outlined in Sect. 3.3 using 1.3-mm and 3-mm 
flux density data only. We compensated the spectral time delay 
by artificially retarding the 1.3-mm lightcurve by 55 days. We 
present the resulting modified spectral index evolution in Fig.fTTI 
(left panel). Time bins are the same as in Fig. |4] (bottom left 
panel) which should be compared to Fig.[TT](left panel). We note 
that retarding the 1 .3-mm data removes the inverted values for a 
that are visible in Fig. |4] (bottom left panel) during the outburst 
2005/2006. 



4.6. 3C 84 

The 3-mm (1.3-mm) lightcurve of this Seyfert 2 galaxy varies 
smoothly during the entire monitoring time between ai3.5 Jy 
(a2 Jy) in 2002 and all Jy (a7 Jy) in 2010. Since about 2004 
the fluxes raise continuously with exponential rise time scales 
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of »6 years (^4 years). The flux distributions show somewhat 
complicated double-peak profiles that can be approximated as 
superpositions of three Gaussians peaking at ^4 Jy, siS.S Jy, and 
~9 Jy (~2 Jy, *4 Jy, and ^6 Jy). 

The spectral index is steep, with a » 0.9 (within errors) until 
about 2005 and then drops smoothly to a » 0.5 after 2007. We 
find the flattest spectral index for the highest flux value; however, 
as this is based on a single data point that occupies a large range 
in flux, this observation does not establish a significant trend. 

The power spectra can be described as simple red-noise spec- 
tra, with powerlaw indices /3 ~ 0.7 for the 3-mm andyS ~ 0.4 for 
the 1.3-mm lightcurve, respectively. 

The discrete correlation function computed from the 1 .3-mm 
and 3-mm lightcurves finds a spectral time delay significantly 
different from zero. The most probable time delay is located at 
T * -270... - 330 days, meaning that the 1.3-mm emission pre- 
ceeds the 3-mm emission by almost one year. We note the high 
absolute value of the maximum correlation, DCF{t - -211 d) - 
1 .04 ± 0.05, which is in agreement with the theoretical maximum 
value of 1. However, as already DCF(0) » 0.7, one might won- 
der how large the difference between t x Q and t ^ -300 actu- 
ally is in absolute terms. We therefore tested the impact of the 
time off'set graphically. For this we normalize the 1.3-mm and 
3-mm lightcurves to zero mean and unity standard deviation, in 
analogy to the computation of the unbinned discrete correlations 
according to Eq.|2] In Fig. [10] we accordingly present the quan- 
tity {Sy - {Sv))/cr for each of the two frequency bands; here Sy 
is the flux density, {S y) denotes the time average of the flux den- 
sity, and cr indicates the standard deviation of the flux density. 
In Fig. [To] we present the normalized data separately for the two 
cases T = days (left panel) and r = -300 days days (right 
panel). Indeed the case r = -300 days shows a better agree- 
ment between the 1.3-mm and 3-mm lightcurves than the case 
T = days. This agrees with the result obtained from the corre- 
lation analysis. However, we note that the lightcurves show clear 
trends that are only partially covered by our observations. This 
can mislead a correlation analysis, meaning that the time delay 
we find has to be taken with care. 

Like for 3C 454.3, we probed the impact of the spectral 
time delay on the calculation of the spectral index. We compen- 
sated the spectral time delay by artificially retarding the 1.3-mm 
lightcurve by 300 days. We present the resulting modified spec- 
tral index evolution in Fig. [H] (right panel). Time bins are the 
same as in Fig. [4] (bottom right panel) which should be com- 
pared to Fig. [TT] (right panel). Comparison of the two diagrams 
indicates that removing the spectral time delay systematically 
shifts a towards steeper values by * 0.1. 



5. Discussion 

5.1. Lightcurves and Flux Distributions 

Our analysis reveals substantial - by factors a!2-8 - flux vari- 
ability on timescales of several years in all sources of our sam- 
ple. Even though all lightcurves (see Fig. [TJ show strong varia- 
tions, the characteristic (exponential rise/decay) timescales dif- 
fer by more than one order of magnitude, ranging from <0.5 
years (3C454.3) to ^7 years (0923-H392). Accordingly, we find 
lightcurves with smooth morphologies during the entire « 14- 
year monitoring time (0923+392, 3C 84) as well as lightcurves 
that are "spiky" (3C 1 11, 3C 454.3). 

The flux distributions (Figs. [2] [3]l of five out of six (the ex- 
ception being 3C 111) AGN show bi- or even multimodality. 



This indicates that the flux variability of the concerned sources 
cannot be described by uniform stochastic processes: 

- In case of Gaussian (white) noise processes the distribution 
of amplitudes would be (obviously) Gaussian. 

- In case of red noise processes (see also Sect. 3.4), the dis- 
tribution of amplitudes is Gaussian for a given sampling fre- 
quency interval [/, / -i- df] (e.g. Timmer & Konig 1995}. 
The distribution of amplitudes resulting when summing over 
all / would thus a priori be a superposition of Gaussians 
of varying widths (i.e. a unimodal, symmetric distribution). 
However, for physical lightcurves one needs to take into ac- 
count that fluxes are constraint to positive values. This intro- 
duces a bias towards higher flux values, meaning that the re- 
sulting flux histogram shows an extended wing towards high 
values. 

- In case of multiplicative stochastic processes the distribution 
of amplitudes would follow a unimodal log-normal distribu- 
tion, i.e. an asymmetric distribution with an extended wing 
towards high values (see Uttley et al. 120051 for a discussion 
of AGN X-ray lightcurves). 

The flux distributions of 3C 1 1 1 are (at least marginally) consis- 
tent with log-normal distributions, pointing towards an underly- 
ing multiplicative process. Such processes can be understood as 
chains of events where each event is triggered by the previous 
one with a certain probability < 1 (e.g., Montroll & Shlesinger 
119821) . Candidates are outbursts that trigger secondary emission, 
or shocks in jets that trigger secondary shocks or fragmentation 
of emission regions. 

For the remaining five sources, it appears to be necessary 
to distinguish separate states of source activity like "quiescent" 
and (maybe multiple) "flare states" where each state is character- 
ized by an individual flux distribution. The final, observed flux 
distribution would then be a superposition of the individual dis- 
tributions. 



5.2. Spectral Indices 

All sources of our sample show (mostly) fairly steep spectral 
indices within ranges a * 0.5 - 1 (Fig. |4]|. All sources show a 
significant, smooth variability of their indices on timescales of 
years. Only 3C 454.3 provides an exception from both of the 
aforegiven statements: it shows a faster variability on timescales 
of »l-2 years and reaches flat and even inverted indices a ^ 
-0.5 - 0. As our sample is composed of very luminous AGN, 
we can safely assume that their mm-emission is dominated by 
synchrotron radiation. We can thus classify the sources of our 
sample (with the possible exception of 3C 454.3) as outflow (i.e. 
halo or jet) dominated sources for which typically a ^ 0.5 - 1 
(e.g. Krolik 1999). This is consistent with the fact that for all 
sources luminous halos or jets on kilo-parsec scales have been 
observed (see Table[TJ. 

The variations of the spectral indices point towards varia- 
tions of intrinsic source properties on timescales of years. A pri- 
ori, there are three potential mechanisms: 

1 . Variations of the optical depth. Synchrotron emission from 
optically thin emission regions shows steep spectral indices 
with a > 0.5; strong absorption in emission regions (i.e. op- 
tically thick emission) leads to approximately flat (a ~ 0) 
or even inverted, complex spectra (see, e.g., Kembhavi & 
Narlikar 1999 for a review). Accordingly, variations of the 
optical depth can lead to variations of the spectral index in- 
between these extreme cases. 
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Table 2. True spectral indices yStme vs. observed indices ySobs- 
Observed indices are given for three cases: no white (Gaussian) 
noise added to the lightcurve (cr - 0), Gaussian noise with cr of 
10% of the flux values added (cr - Q.lSv), Gaussian noise with 
cr of 16% of the flux values added (cr = Q.l6Sv)- These values 
are obtained from Monte Carlo simulations. Errors indicate the 
Gaussian widths of the distributions. 



PtiLie 


ySob.(o- = 0) 


yS„b,(o- = 0.15.,) 


P,U(r = 0A6S,) 


2.0 


1.84±0.11 


1.46±0.23 




1.33±0.25 


1.5 


1.51±0.16 


1.30±0.14 




1.17+0.17 


1.0 


0.99±0.13 


0.95±0.11 




0.85±0.13 


0.5 


0.48±0.15 


0.44±0.13 




0.39±0.10 



2. Variations of the relative luminosities of different source 
components. Emission from (optically thick) nuclei shows 
usually fairly flat spectra {a ^ 0) whereas the emission from 
(optically thin) outflows (jets, lobes) typically shows steep 
spectra (a « 0.5 - 1; e.g. Krolik [19991 ). Variations in the rel- 
ative luminosities of the various components therefore mod- 
ify the spectral index of the source-integrated emission. The 
spectrum becomes steeper if the relative contribution by out- 
flows increases, it becomes flatter if the relative contribution 
by the nucleus increases. 

3. For optically thin synchrotron emission from an ensemble 
of electrons with a powerlaw energy spectrum with index 
y the spectral index of the radiation is given by a = (y - 
l)/2 (e.g. Ginzburg & Syrovatskii 1965); for a = 0.5 - 1 
this corresponds to y = 2 - 3. Accordingly, variations in the 
electron energy spectrum express themselves as variations of 
the spectral index. 

For three of our six target sources (namely 3C 1 11, 3C 273, 
and 3C 454.3) the spectral indices are - at least occasionally - 
significantly flatter than the range permitted for the assumption 
of optically thin emission (i.e. we observe a < 0.5) for parts 
of the monitoring time. This suggests that scenario 3 cannot be 
applied regularly as this scenario is valid only for optically thin 
emission. 

Additional information is provided by the relation between 
source fluxes and their spectral indices (Fig. |5]). Especially, the 
Sy - a relation reveals that scenario 2 can play only a minor 
role for the variability of the spectral index: if the a variations 
were due to (or at least dominated by) luminosity increases or 
decreases of parts of the sources, there should be a tight corre- 
lation between source fluxes and spectral indices - which we do 
not observe. From the combined information we therefore con- 
clude that variations of the optical depth with time (scenario 1) 
are the most likely mechanism for variations of the spectral in- 
dex. 



5.3. Power Spectra 

The power spectra of our sources (Fig. |6] |7]i globally follow 
red noise distributions, meaning powerlaws with spectral indices 
jS > (cf. Section 4.3). This, as well as the apparent absence of 
periodic signals, shows that the emission processes are intrin- 
sically stochastic. The power-spectral indices /3 « 0.4 - 1 are 
located between those of Gaussian (white) noise (/3 = 0) and 
flicker noise (J3 - 1) processes (e.g. Timmer & Konig [T995I I. 
They are significantly flatter than the indices found from optical 
and cm/radio (e.g. Hufnagel & Bregman 1992 ) as well as X- 
ray (e.g. Lawrence & Papadakis |1993l l studies of AGN emission 
which are located in the range 1-2. 



When discussing the power spectral indices we have to take 
into account that we observe lightcurves that are affected by 
white (Gaussian) measurement noise. This noise is dominated 
by a systematic scatter that amounts to a; 10% for the 2-mm and 
3-mm fluxes and ^16% for the 1.3-mm flux data (see Sect. 2). 
The presence of white measurement noise leads to the obser- 
vation of power-spectral indices /3 that are actually flatter (i.e. 
closer to the case of white noise, /3 - 0) than the indices intrinsic 
to the source emission. In order to estimate the intrinsic spectral 
indices, we performed Monte Carlo tests using simulated red- 
noise lightcurves with known ideal spectral indices. For this test 
we applied the algorithm by Timmer & Konig ( 119951 ). We calcu- 
lated artificial lightcurves coiTesponding to red noise lightcurves 
with /? = 2, y6 = 1.5, yS = 1, and y8 = 0.5, respectively. For 
each case, we examined three scenarios: the absence of mea- 
surement noise, Gaussian measurement noise with cr = Q.lSv 
(i.e. a Gaussian width corresponding to 10% of the observed 
fluxes), and Gaussian measurement noise with cr - 0.165v (i.e. 
a Gaussian width coiTesponding to 16% of the observed fluxes). 
We simulated 100 lightcurves with 500 data points each for each 
of the 12 cases examined. In order to assess the impact of gaps 
in the data, we applied a sampling approximately correspond- 
ing to the worst case actually encountered which is the sampling 
of the 1.3-mm lightcurve of 3C 111 (compare Fig. [1] top right 
panel). For this, we inserted - by simply leaving out data - a 
total of eleven artificial gaps, each covering ^2% to ^9% of the 
total observing time. For each artificial lightcurve we calculated 
a Scargle periodogram (Eq. |6]l and derived the powerlaw index. 
We present the results in Table|2] Due to the stochastic nature of 
the noise signals, the "observed" values for /? show substantial 
scatter around the averages; the Icr widths of the coiTesponding 
distributions are quoted as uncertainties in Table |2] 

From our test we can draw various conclusions: 

- The presence of Gaussian measurement noise flattens the 
power-spectral indices to^Sobs < 1.5 even ifySn-ue ~ 2. 

- AsySobs ~ 0.4-0.7 for all sources except 3C 111, the intrinsic 
power-spectral indices should be flatter than 1 in general. 

- Increasing the amplitude of measurement noise from Q.lSv 
(like for our 3-mm flux data) to 0.165'v (like for our 1.3-mm 
flux data) leads to a further flattening of the power-spectral 
index by about 0.1. This agrees with our observations: we 
find a systematic trend towards flatter (by ^0.1) indices for 
the 1 .3-mm data compared to the 3-mm power spectra of the 
corresponding sources. 

This also shows that the diff'erence between our results and 
those reported by Hufnagel & Bregman (1992) or Uttley et 
al. (2002) cannot be attributed to measurement noise effects. 
Instead, we have to attribute this difference to the different fre- 
quency regime (and thus emission region) we observe. This in- 
dicates that at millimeter wavelengths, adjacent AGN emission 
events tend to be less correlated than in other regimes. 

The fact that intrinsically y6 < 1 (except for 3C 111) points 
towards stochastic emission processes located between the cases 
of white noise and flicker noise, with 3C 1 1 1 being closer to the 
case of random walk noise (/? = 2). Formally, the various cases 
are connected by different levels of correlation among adjacent 
data. In case of white noise (/? = 0) adjacent data are completely 
uncorrelated. Random walk noise (J3 - 2) corresponds to the 
integral of white noise over an infinite interval in time, mean- 
ing that adjacent data are fully correlated. Intermediate cases 
(0 < yS < 2) can be expressed as "half" integrals of white noise 
(again, over infinite intervals in time). Such a "half" integration 
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is performed by convolving a white noise time series with a pow- 
erlaw kernel G oc r^ with ^ > and t being the time. For the 
specific case of flicker noise (/?=!) this was discussed in detail 
by Press (I1978I I. leading to ^ = 1 /2. From this, one may a pri- 
ori expect to find any < ^0 < 2 in AGN time series driven by 
stochastic emission processed 

Even though there is a straightforward mathematical recipe 
for the creation of red noise (in the range < /? < 2) time se- 
ries, there is not yet a consistent understanding which physical 
mechanisms create those time series. This is a problem of im- 
portance not only in astrophysics but in many fields of physics, 
and various approaches towards a general explanation have been 
proposed (e.g., Press ll978l Pang et al. ,1995. M akse et al. ll996l l. 
More recently, Kaulakys & Alaburda ( I2009I I and Kelly et al. 
(1201 1) proposed time series models based on superpositions 
of exponentially decaying outbursts occuring at random times 
and with random amplitudes. Both studies are able to reproduce 
the occurence of powerlaw power spectra. However, as already 
noted by Press d 19781 1 who discussed (and eventually rejected) 
similar proposals, those approaches might only shift the prob- 
lem from the question "What process creates those power spec- 
tra?" to the question "What process generates those outburst se- 
quences?" without actually solving the problem. 

Although the power spectra globally follow red noise laws, 
tiiree out of six sources (3C 111, 3C 273, and 3C 454.3) ai-e 
hardly consistent with a single stochastic emission process. 
Instead, we observe excess (by factors a!2-3) spectral power 
with respect to the average powerlaw for sampling frequen- 
cies / < 0.005 - 0.01 days '. In combination with the in- 
formation from the flux distributions (Sect. 5.1; Fig. |2]|3]l this 
strongly suggests that "quiescent" high-frequency and "flaring" 
low-frequency emission have to be treated as separate source 
states. This conclusion deviates from the more simple picture 
drawn by other studies (e.g., Hufnagel & Bregman 1992; Uttley 
et al. I2OO2I 1 that treat AGN time series and power spectra as 
uniform red-noise distributions. However, a physical dichotomy 
of quiescent vs. flaring source states has been suggested for the 
near-infrared emission observed from the Galactic nuclear black 
hole, Sagittarius A* (e.g. Gillessen et al. 120061 Trippe et al. 
120071 Dodds-Eden et al. 1201 11 1, and has been hotly debated since 
then (compare, e.g., Meyer et al.l 2007l vs. Meyer et al.'2008). 

We also note that a segregation of distinct source states is 
well-known for the case of Galactic microquasars. These sources 
mirror many properties of AGN (taking into account scaling 
with the black hole mass) and their emission shows distinct 
states that are ususally classified as "hard" and "soft" accord- 
ing to their X-ray spectral properties (e.g., Tananbaum et al. 
WnX Mirabel & Rodriguez 19981 [19991 Takeuchi & Mineshige 
11998 ). A possible explanation for this transition (Klein- Wolt et 
al. |2002) is the difference between the ejection of fairly localized 
knots from the black hole vicinity (leading to X-ray and radio 
flares) and the ejection of quasi-continuous, extended jets that 
are initially optically thick (leading to "plateaus" in the X-ray 
and radio lightcurves). For micro-quasars, these events happen 
on timescales from minutes (minimum time between the ejection 
of knots) to tens of days (maximum length of "plateau" phases 
observed); for AGN, corresponding phenomena would have to 
occur on timescales of years - thus deviating strongly from a 



simple Unear black-hole mass scaling relation - in order to be 
consistent with our observations. 

Notably, the "flare" and "high-frequency" regimes of the 
power spectra are consistent with having the same powerlaw in- 
dex though on different (by factors a!2-3 in spectral power) en- 
ergy levels. This indicates that both emission regimes, although 
they appear to be physically disconnected, are caused by intrin- 
sically stochastic processes with similar correlation properties. 
Remarkably, the three sources in question - 3C 111, 3C 273, 
and 3C 454.3 - are also those for which we find the (temporar- 
ily) flattest spectral indices, indicating optically thick emission. 

5.4. Spectral Time Delays 

We tested the presence or absence of time lags between the 1.3- 
mm and 3 -mm fluxes by means of a correlation analysis (cf. 
Eq. [8]l. We can exclude the presence of significant - referring 
here to a 3cr confidence level - offsets in time for four of our 
tai-gets, namely 0923-H392, 3C 111, 3C 273, and 3C 345 (see 
Fig-Ejl. This statement holds within accuracies (referring to the 
ranges of timelags consistent with a zero offset within the 3cr 
confidence band) of few tens of days. Our results are consistent 
with those obtained by the time delay analyses of Tornikoski 
et al. (1 19941 1 and Hovatta et al. (2008) in the cm-radio range. 
Hovatta et al. (2008) conclude that those limits on time delays, 
both in the cm- and mm-regimes, can be understood in the frame 
of the semi-empiric "generalized shock model" by Valtaoja et al. 
(119921) . 

For the remaining two sources, 3C 454.3 (after 2005.2; cf. 
Fig.|9]) and 3C 84, we do find indication that the 1.3-mm emis- 
sion preceeds the 3-mm emission by » 55 days (3C 454.3) and 
~ 300 days (3C 84). We note that the result for 3C 454.3 is con- 
sistent with the results by Tornikoski et al. ( 1994) and Hovatta et 
al. (2008 ) who do not find a time delay in their radio lightcurves; 
the data of Hovatta et al. ( 120081) end in 2005.3, just at the begin- 
ning of the "flare phase" of 3C 454.3. 

For 3C 454.3, we suspect variations of the optical depth with 
frequency as the principal cause for the time delay. Adiabatically 
expanding plasmas become optically thin at higher energies first, 
introducing a systematic time delay between two lightcurves ob- 
tained at two frequencies with the high-frequency emission lead- 
ing. Details depend on various parameters, notably the electron 
energy distribution, the optical depth at a reference frequency, 
the magnetic field strength, the electron density, and the emis- 
sion region size (van der Laan |1966l ). This mechanism has been 
incorporated into the shock-in-jet model by Marscher & Gear 
( 1985 ); proper quantitative treatment takes into account the com- 
bination of Compton, synchrotron, and adiabatic cooling (e.g., 
Fromm et al. 2011', and references therein). Our interpretation 
is supported by the fact that the spectral index shows fast vari- 
ability from inverted to almost optically thin spectral indices in 
2005-2007 and remains in the optically thick regime after 200'cl 
This is in agreement with the results by Jorstad et al. (120101) who 
find time delays up to a!215 days between 37 GHz and 230 GHz 
lightcurves, with the 230 GHz emission leading. 

For 3C 84, the answer is less clear High-resolution radio- 
interferometric observations let Asada et al. 6200 6) conclude that 
the lightcurve, as observed until 2001, is in agreement with adi- 
abatic cooling of an expanding radio lobe. The rise of the emis- 



* Altiiougii there seems to be some confusion in the literature; e.g., 
Do et al. ( 2009 1 conclude that the infrared emission from Sagittarius A* 
follows a red noise power spectrum with /3 ~ 2.5, without providing a 
physical interpretation of this unexpected result. 



'' This observation changes only marginally when artificially remov- 
ing the spectral time delay; see Fig. \TT\ The most notable difference 
between the two calculations is the time when inverted indices are ob- 
served. 
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sion around 2005 nicely agrees with observations of the ejection 
of a new radio jet component in 2005 (Nagai et al. l2010l l. Those 
observations are a priori in agreement with optical depth varia- 
tions in the context of the shock-in-jet model. However, we note 
the unusually large time delay of a;300 days - assuming the re- 
sult of the DCF analysis is correct. Such a large delay suggests 
additional travel time effects given by the geometry of the emis- 
sion region. If the emission is composed of various components 
emitted from spatially separate zones, according time lags may 
occur. Given the observed time offsets of siSOO days, this would 
suggest a spatial extension of the emission region of the order of 
one light-year. 

6. Conclusions 

We have analyzed the long-term activity of six luminous active 
galactic nuclei using mm/radio lightcurves spanning ^14 years 
in time. From our analysis we draw the following principal con- 
clusions: 

1. All sources show substantial, non-periodic variabiUty by fac- 
tors ^2-8 on timescales of years. This is in good agreement 
with the results reported by Hovatta et al. ( 120071 IJ50 8 ) and 
shows that source monitoring times on scales of decades 
are necessary in order to obtain a realistic overview on the 
source activity. We also note that at least for the special 
case of Sagittarius A* both observational (e.g. Koyama et al. 
120081) and theoretical (e.g. Cuadra et al. l2008l l evidence sug- 
gests activity timescales of several hundred years. An even 
more extreme example is provided by recent observations 
of "Hanny's Voorwerp" which suggest flux variations of the 
galactic nucleus of IC 2497 by factors >100 on timescales 
<70,000 years (e.g., Schawinski et al. 2010). Therefore one 
should be aware that AGN monitoring studies, like our work, 
might observe some of their targets at special times even 
when covering timelines of decades. 

2. The flux density distributions of five out of six sources show 
clear signatures of bi- or even multimodality. This indicates 
that the emission processes cannot be treated as uniform 
stochastic processes. Instead, we have to assume a segre- 
gation of "quiescent" and (maybe multiple) "flare" source 
states that correspond to different physical regimes. Each 
state would be characterized by an individual (probably red- 
noise type) flux distribution. 

3. Our sources show mostly steep (a ^ 0.5 - 1) spectral indices 
that indicate outflow (halos, jets) dominated synchrotron 
emission. The spectral indices show substantial variability 
that is not correlated with the source fluxes. From the com- 
bined information we conclude that optical depth variations 
are the dominant mechanism for the spectral variability. 

4. All sources show power spectra that are globally consistent 
with red noise (i.e. /? > 0). We find that the intrinsic (i.e. mea- 
surement noise corrected) power-spectral indices are flatter 
than 1 for five of our sources, placing their emission between 
the cases of flicker noise (J3 - I) and white noise (J3 = 0). 
The only exception is 3C 1 1 1 whose power spectrum is con- 
sistent with being located between flicker noise and random 
walk noise (J5 - 2). For three sources we find that the low- 
frequency ends of their power spectra appear to be upscaled 
in spectral power by factors =s2-3 with respect to the overall 
powerlaws, presenting another indication towards a segrega- 
tion between "quiescent" and "flaring" source states. 

5. For four of our targets - 0923H-392, 3C 111, 3C 273, and 
3C 345 - we find that the 1.3-mm and 3-mm lightcurves are 



simultaneous within 3cr confidence levels. This is in agree- 
ment with shock-models like the "generalized" model by 
Valtaoja et al. (1992); given the measurement accuracies of 
tens of days however these constraints are fairly loose. For 
3C 454.3 (in its "flai-e state(s)" after 2005.2) and 3C 84 we 
find spectral time delays of s!55 days and s; 300 days, respec- 
tively, with the 1.3-mm emission leading. The behaviour of 
3C 454.3 can be traced back to variations of the optical depth 
with time and frequency, in agreement with previous work. 
In case of 3C 84, we suspect a combination of optical depth 
and emission region geometry effects as cause for the time 
lag. 
6. We do not see a correlation between the observed source 
properties - variability, spectral index, power spectrum - 
with physical source parameters - redshift, black hole mass, 
kilo-parsec morphology - as given in Table [T] 

The discovery of distinct source emission states challenges 
the common assumption that AGN emission can be described by 
uniform stochastic processes. A promising road for further stud- 
ies is indicated by the cases of microquasars, Sagittarius A*, and 
IC 2497 where similar segregations of emission states have been 
seen. Long-term, multi-wavelength monitoring of AGN appears 
to be the most important tool towards a deeper understanding of 
this phenomenon. 
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Fig. 1. The long-term lightcurves of our six target quasars spanning timelines from 1996 to 2010. Please note the different flux axis 
scales. For each source we give the fluxes vs. time separatly for the 1.3-mm (here labeled 1-mm for brevity), 2-mm, and 3-mm 
bands. For each source (except 3C 111) we approximate the most prominent, continuously sampled variations as exponential rises 
or decays according to 5v oc exp(±f/ff) (orange lines). Here t is the time, fj the characteristic rise/decay timescale of the curve 
labeled x (x - a, b, c, ...). The sign of t is +t for exponential rises and -t for decays. All parameters are given in the corresponding 
diagrams. 
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Fig. 2. Distributions of the 3-mm (left column) and 1.3-mmband (right column) fluxes of 0923+392, 3C 111, and 3C 273. Please 
note the different axis scales and bin sizes. In order to minimize the impact of irregular sampling, we binned all lightcurves (compare 
Fig. [1]) in time using A^ bins with sizes Af. The parameters A^, Af, and the sizes of the flux bins ASv are given in the respective 
diagrams. Errorbars indicate binomial errors. Continuous grey lines indicate approximations to the data composed of individual 
Gaussians (dashed grey lines). In three cases (0923+392 at 1.3 mm, 3C 111 at 1.3 mm and 3 mm) we also indicate the best-fitting 
log-normal distributions (dotted black curves). 
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Fig. 3. Same as Fig.|2]for 3C 345, 3C 454.3, and 3C 84. 
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Fig. 4. Evolution of the spectral indices from 1996 to 2010. Please note the different axis scales. The spectral index a is defined via 
Sy cc v^" . Errorbars along the time axes denote the bin sizes, error bars along the a axes denote the statistical Icr errors. Horizontal 
dashed lines (where visible) mark the a = lines. For each source we chose the bins such that each bin covers a phase of similar 
flux levels (quiescent phases, flares). 
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Fig. 5. Evolution of the spectral index a as function of the 3-mm flux. The data points are connected by dotted lines according to 
their order in time. Labels 1, 2, 3, ... mark the first, second, third, ... measurement in time. The time bins used here are those defined 
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Fig. 6. Scargle periodograms (spectral power Af vs. frequency /) of the 3-mm (left column) and 1.3-mm band (right column) 
lightcurves of 0923+392, 3C 111, and 3C 273. Please note the different axis scales. The black curves denote the measured power 
spectra, the continuous grey lines marks the best-fitting powerlaws. The spectral index /? of the power spectra is defined via Aj- oc /''. 
Dashed grey lines (where present) correspond to the best-fiting powerlaws (continuous grey lines) upscaled by factors s as given in 
the corresponding plots. 
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Fig. 7. Same as Fig.|6]for 3C 345, 3C 454.3, and 3C 84. 
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Fig. 8. Discrete correlation functions (DCF) as functions of time lags t between the 1.3-mm and 3-mm lightcurves (see Fig.[lJ of 
each AGN. The black graphs indicate the measured DCF, bin sizes At are indicated in the corresponding diagrams. The maxima 
of the curves correspond to the most probable time lag between the spectral bands. Grey shaded areas mark the DCF ranges 
[max(DCF), maxiDCF) - 3cr] with cr being the statistical error of max(DCF). Thus all values located within the shaded area are 
consistent with the maximum of the DCF curve within 3cr confidence. A positive (negative) time lag means that the 3-mm lightcurve 
preceeds (lags behind) the 1.3-mm lightcurve. Vertical dashed line mark the null positions of the time lag axes, corresponding to the 
absence of time lags. 
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Fig. 9. Like Fig. [8] for 3C 454.3 before 2005.2 (left panel) and after 2005.2 (right panel). During the "quiescent" phase before 
2005.2, the peak of the DCF is in good agreement with a time delay of zero. We note the low absolute value of the correlation, 
DCF(Q) X 0.6. For the "flaring" phase after 2005.2, we find a significant time delay in the range (as defined by the 3cr confidence 

interval) of 15... - 80 days, peaking at r a; -55 days. The negative delay indicates that the 1.3-mm lightcurve is ahead in time 

of the 3 -mm lightcurve. 
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Fig. 10. A visualization of the time delay between the 1.3-mm and 3-mm lightcurves of 3C 84. Left panel: The 1.3-mm and 3-mm 
lightcurves. Each lightcurve is normalized to zero mean and unity standard deviation, in analogy to the calculation of the unbinned 
discrete coiTelations according to Eq. Q Here S y indicates the flux density, {S y) is the average flux density, and cr denotes the 
standard deviation of each lightcurve. Right panel: Same as left panel, but with the 1.3-mm lightcurve artificially retarded by 300 
days. Comparison of the two panel shows improved agreement of the two lightcurves when taking into account the time delay 
T X! -300 days found from the discrete correlation function; see also Fig. |7l bottom right panel. 
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3C454.3, 1.3— mm lightcurve retarded by 55 days 3C84, 1.3— mm lightcnrve retarded by 300 days 



Fig. 11. Spectral index a vs. time based on 1.3-mm and 3-mm flux data only, with the U-mm lightcurves artificially retarded. Left 
panel: 3C 454.3, using a spectral time delay t = -55 days. Right panel: 3C 84, using a spectral delay r - -300 days. The time 
delays are given by the correlation analysis illustrated in Fig.|7] This figure should be compared to Fig.|4] 



20 



